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CO ' Abstract The drastic rise of the proton structure function F 2 (x, Q ) 

when the Bjorken variable x decreases, seen at HERA for a large span 

q ! of Q 2 , negative values for the 4-momentum transfer, may be damped 

when Q 2 increases beyond ~ several hundreds GeV 2 . A new data anal- 
ysis and a comparison with recent models for the proton structure func- 
tion is proposed to discuss this phenomenon in terms of the derivative 



g* d£nF 2 (x,Q 2 )/d£n(l/x) 



1 Introduction 



The so-called "HERA effect" discovered some years ago by the experimentalists in deep 
inelastic scattering (DIS) and anticipated by a number of theoreticians has arisen a great 
interest in the community (see e.g.]!] and references therein). It concerns the strong rise of 
the proton structure functions (SF) F 2 (x, Q 2 ) when the Bjorken variable x decreases, for the 
experimentally investigated Q 2 , negative values for the 4-momentum transfer. 

A slow-down of a further rise is inevitable, as a consequence of unitarity, if however the 
Froissart-Martin 0] bound is valid for j*p interaction at least at x — > 0. In this case a power- 
like behavior of SF at small-x must be transformed into a logarithmic-like one (like in hadron 
amplitudes, when a procedure of unitarization or eikonalization is applied). Evidently, such 
an evolution leads to a damping of the fast growth of SF at x — > 0. 

A natural question arises : will this HERA effect still subsist in any x — Q 2 region to be 
investigated, or more pragmatically in which kinematical region a damping is to be expected ? 

It is not easy to exhibit this interesting tendency in the behavior of F 2 (x,Q 2 ) at small 
x and rising Q 2 . This effect is very weak because a sufficient amount of data at high Q 2 is 
still lacking. Furthermore, it may stand near the kinematical limit at the available energies 

1 E-mail: desgrolard@ipnl.in2p3.fr 
2 E-mail: sasha@len.uzhgorod.ua 
3 E-mail: martynov@bitp.kiev.ua 



(it is the case at HERA [|, thus preventing a clear evidence on the basis of existing data. 
Nevertheless, one can guess it by simply inspecting by eye the evolution of the icons showing 
the experimental SF versus x at given Q 2 increasing up to the highest values (see in Fig.l 
such a representation using a double logarithmic scale); if this inspection seems inconclusive, 
one can at least insure that no contradiction appears. 
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Figure 1: Experimental data of ZEUS and HI collaborations at high Q^ 



Aside from this empirical mean, a way to quantify the preceding effect is to consider the 
derivative of the logarithm of the proton SF with respect to £n(l/x) (x— slope for brevity) 
defined as 



B x (x,Q 



, d£nF 2 (x,Q 2 ) 



(1) 



d£n(l/x) ' 

the damping of the HERA effect would correspond to the presence of a maximum when 
plotting B x versus Q 2 for a given low x. 

It is worth noting that this derivative is related to an "effective power" (exponent of 1/x) 
A, sometimes labelled A, introduced when redefining the proton SF as 



F 2 (x,Q 2 ) = G(Q 2 



nA(x,Q 2 ) 



X 



(2) 



4 thc kinematical limit y = —r^- — tt < 1 for HERA measurements, with s — mi ~ AE e E„, in terms of the 

y x(s—nrip) — ' P e y' 

positron E e and proton E p beam energies, writes Q 2 ( GeV 2 ) < 10 5 a; since August 1998. 



or in other terms to a Pomeron "effective intercept" oup (depending on x and Q 2 ) 

a e v ff (x,Q 2 ) = l + A(x,Q 2 ). (3) 

We emphasize that we can identify B x = A, only in the case of a model where A is x- 
independent; in general B x and A do not coincide. A careful analysis of experimental data, 
based on a comparison of various parametrizations of the SF, has suggested [[J a possible 
slowing down of the SF rise revealed as a growth of F 2 with decreasing x becoming steeper 
and steeper (when Q 2 rises) which is changed at Q 2 ~ 200 — 500 GeV 2 into a growth 
saturating and even becoming less and less steep at further increasing Q 2 . 

From a phenomenological point of view, many parametrizations of the proton structure 
functions successfully accommodated for the whole - or a part of - available data set and in 
particular for the steep rise of the SF when x decreases for a large span of Q 2 values (see 
e.g. [4-17]). Several of them have given explicitly some hints (see e.g. [[[0], [16], |T7|) on the 
existence of a slowing down of the HERA effect. 

We believe that the region of the contour x — Q 2 plane with x < 10 _1 and Q 2 > 100 
GeV 2 is of particular interest for the above question. Happily, recent measurements of SF 
have become available, from HI and ZEUS collaborations which fulfil partly this kinematical 



criterion. They complete or correct the previous data near the HERA collider |Tj| |19|, g0[ ^T 
and P, G3L M E3] and from fixed target experiments EB, 1271, ER ESI BU 



Here, we actualize the conclusions of previous paths driving towards the existence of a 
damping effect. We perform with new data a relevant analysis similar to |3j, where, however, 
only the averaged slopes B x (x,Q 2 ) were determined in almost the whole available interval 
of x at given fixed Q 2 . 

It is necessary to remark here that the data on A e // [EJJ cannot be interpreted as showing 



a dependence of the effective intercept (A e jf(Q 2 )) on Q 2 . Each point B x has not been 
extracted at the same < x >, There is an evident correlation between the values of Q 2 and 
of the corresponding < x > , namely a larger < x > corresponds to a larger Q 2 for almost all 
points f\. To exhibit a true dependence of the effective intercept on Q 2 , one need to have the 
local B x extracted at the same < x > but for different Q 2 . It is one of the reasons motivating 
the present work. We attempt a new, more detailed, systematic (on all available data) and 
model independent analysis of the local slope B x (x, Q 2 ) having in mind its possible variation 
with x. It allows not only to extract the averaged slopes at fixed Q 2 but also to determine 
the dependence of B x (x, Q 2 ) on Q 2 at fixed x. 

We also turn towards three basically different phenomenological models for the proton SF, 
namely: the "ALLM" model ||, the "interpolating Lyon-Kiev-Padova" (LKP) model [|T^], 
the "soft dipole Pomeron" (DP) model |17j, and have been refitted (or eventually improved 



as concerns LKP) with a common set of data including the most recent ones, their predictions 
concerning the x— slope are compared to the results of the present analysis. 

2 Models for the proton structure functions 

2.1 Common set of data included in the fitting procedure 

The choice of the data set used in the fits may have crucial consequences when calculating the 
derivatives of the SF, especially when this quantity is to be discussed near the experimental 
kinematical limit or/and near the theoretical validity limit. For these reasons, a unique set 

5 See also a discussion about A e ff in the second paper of M. 



including the most recent data has been used in new fits of three models of the proton SF 
we choose to discuss. 

These updated data are listed and referenced in Table 1. The complete available experi- 
mental kinematical range is taken into account for Q 2 (i.e. 0. < Q 2 < 30000 GeV 2 ) and x 
(i.e. 6. 10~ 7 < x < 0.75), with the lower limit of the center of mass energy for 7*p scattering 
W > 3 GeV. 



Observable 
Exp. — year of pub., 



Rcf 



Nb points 



X 
model [16] 



X 
model [17] 



Fl 

HI - 1995 

HI - 1996 

HI - 1997 

HI - 1999 

ZEUS - 1996 

ZEUS - 1997 

ZEUS - 1999 

ZEUS - 2000 

NMC - 1997 

E665 - 1996 

SLAC - 1990/92 

BCDMS - 1989 



93 
129 
44 
130 
188 
34 
44 
70 
156 
91 
136 
175 



76.1 
90.9 
53.0 
237.9 
259.9 
12.2 
26.7 
139.6 
285.5 
109.7 
193.6 
274.8 



70.7 
68.9 
53.5 
111.0 
261.5 
17.3 
38.5 
79.3 
161.1 
99.5 
107.3 
269.9 



(T 



7.P 
tot 



1975/78; ZEUS - 1994; HI - 1995 |q] 



99 



148.6 



122.7 



Total 



1389 



1907.6 



1461.2 



X 1 I d.o.f. 



1.38 



1.07 



Table 1. Observables sets used in the fitting procedure. Also shown are the distributions of 
the partial x 2 f° r each subset of data used in our new fits with the parameters listed below 



in Tables 2-3, for the two updated models [16, 17 



Generally speaking, we must note a good agreement between both HERA experiments, 
HI and ZEUS, and the fixed target ones. Following the suggestion from |nj, some data from 
[|3~9f| are considered as obsolete and superseded. These correspond to (Q 2 > 250 GeV 2 , for 
all x), (Q 2 = 200 GeV 2 , for x < 0.1) and (Q 2 = 150 GeV 2 , for x < 0.01). 



2.2 Parametrizations of the structure function 

2.2.1 The ALLM model 

The so-called ALLM parametrization [[5] for describing the total 7*p cross section, connected 
to the proton SF, has been updated in 1997 to account for the set of available data in the 
whole range studied so far: < Q 2 < 5000 GeV 2 , and 3.10~ 6 < x < 0.85. Here, to make 
the comparison with the other results more meaningful, we extend the x, Q 2 ranges over all 
the data of Table 1 and perform a new fit of the 23 parameters of this model. We find again 
quite a good agreement which can be characterized by xi o .f=l-06 



2.2.2 The Lyon-Kiev-Padova model revisited 



We emphasize that this LKP model was created for interpolating between soft and hard 
dynamics, combining Regge behavior and the high Q 2 asymptotics of the DGLAP evolution 
equation, for the low— x range, with a small number of free parameters (8, in fll6|). 

To improve the quality of the fit on all data, we accept to lose the simplicity of the 
original model by modifying here the " large" x extension which is no more borrowed from the 
relatively simple and rather efficient model M but, instead, we follow the more sophisticated 
approach from W%. Explicitly, the SF has been rewritten 



F 2 (x,Q 2 ) = Ff < 0) (*,Q 2 ) • (1 - x) p ^ + Ff 5 < 0) (*,Q 2 ) • (1 - x) R ^ 



(NS,0) . 



(4) 



where the Pomeron contribution (singlet Fl \ ) and the effective secondary Reggeon com- 
ponent (non-singlet _F 2 ) have here the same expressions as in |n| at low— a;, and with 
the following Q 2 -dependent exponents of the large— a: factors 



P(Q 2 



Poo + 



PO -Poo 

i + QVQl 



R{Q 2 



+ 



r 



1 + QVQ 2 r 



(5) 



where p^, po, r^, ro, Q p , Q r are thus additional parameters in the present version of the LKP 
model . In definitive, taking into account the large x in practice doubles the parameters 
number, but it is the price of such an extension. The details on the performances of this 
extended model are given in Table 1 in terms of the \ 2 f° r each subset of data, yielding 
Xd *=1.38. The parameters we find are listed in the Table 2, with their meaning explained 
in [lq| and above. 



Parameters 


previous work [16] 


present work 


low x— terms 






A 


0.1470 


0.1190 


a (GeV 2 ) 


0.2607 


0.2300 


e 


0.08 (fixed) 


0.0895 


72 


0.0200 


0.0221 


Ql (GeV 2 ) 


0.1675 


0.1946 


Q\ (GeV 2 ) 


1174. 


7800. 


B 


0.7575 


1.6409 


b (GeV 2 ) 


0.4278 


1.46 


a r 


0.5241 


0.48 (fixed) 


large x— terms 


(fixed) 




Q' 2 P (GeV^) 


- 


1.1180 


Po 


— 


«0. 


Poo 


— 


15.093 


Q 2 . (GeV 2 ) 


- 


12.563 


^0 


— 


2.394 


^oo 


- 


3.728 



Table 2. Parameters used in the previous and present versions of the LKP model. The 
non-fitted values are: 71 = 2.4 (suggested by QCD) for both versions, e = 0.08 (from |3"5fl ) 
and the large x exponents (from ||) for the previous work only. 

When the intercept a r of the /-Reggeon is allowed to be free we obtained a better x 2 (~ 
1.22). However the resulting value of a r ~ 0.24 is totally unacceptable from a physical point 
of view. For that reason, we fixed it, choosing a r = 0.48, a "mean" value used in the models 
with a supercritical Pomeron. 



2.2.3 Soft dipole Pomeron model 



The refitted parameters of this model are given in Table 3 together with the old ones pub- 
lished in |T7| (the same notations are used). The most interesting property of this model 
developed in |L7| is the intercept of Pomeron. It does not depend on Q 2 and moreover it is 
equal to one. Nevertheless due to an interference of two Pomeron components, the model 
describes well all observed properties of the proton structure function, including a growth 
with Q 2 of the effective Pomeron intercept (see for details |]I7[ and [31]). The first Pomeron 
component, dominating at small x, leads to F2(x, Q 2 ) oc £n(l/x) at x — > while the second 
one gives asymptotically a constant and negative contribution to F 2 playing an important 
role in the whole kinematic region where the data on the proton SF data are available. 



Parameters 


previous work [17] 


present refit 


Pi -term 






A* 


.10000E+01 (fixed) 


.10000E+01 (fixed) 


a v {0) 


.10000E+01 (fixed) 


.10000E+01 (fixed) 


9\ (mb) 


.21898E-01 


.22198E-01 


Q\ (GeV 2 ) 


.15400E+02 


.70711E+01 


Q\ d (GeV 2 ) 


.17852E+01 


.14774E+01 


Q\ h (GeV 2 ) 


.33435E+01 


.67975E+01 


dloo 


.13301E+01 


.12601E+01 


dio 


.14370E+02 


.66975E+01 


^loo 


.21804E+01 


.28712E+01 


&10 


.42596E+01 


-.20279E+01 (fixed) 


P2-term 






92 (mb) 


-.99050E-01 


-.10176E+00 


Q 2 (GeV 2 ) 


.34002E+02 


.13748E+02 


Ql d (GeV 2 ) 


.12327E+01 


.15954E+01 


Q 2 b (GeV 2 ) 


.20702E-01 


.80605E+01 


cfooo — dloo 


.00000E+00 (fixed) 


.00000E+00 (fixed) 


^20 


.22607E+02 


.64794E+01 


^2oo 


.24686E+01 


.34510E+01 


^20 


.17023E+03 


.12922E+01 


P-term 






«/(0) 


.80400E+00 (fixed) 


.80400E+00 (fixed) 


9f (mb) 


.29065E+00 


.29405E+00 


Q) (GeV 2 ) 
Q% (GeV 2 ) 


.29044E+02 


.10182E+02 


.54462E+00 


.70413E+00 


Q? b (GeV 2 ) 


.20656E+01 


.84803E+00 


dfoo 


.13554E+01 


.13149E+01 


dfo 


.75127E+02 


.19746E+02 


fc/oo 


.27239E+01 


.33642E+01 


&/0 


.64713E+00 


-.27968E+01 (fixed) 



Table 3. Parameters refitted in the soft Dipole Pomeron model of ||17|| , the original values 
are also quoted. 

Analyzing the properties of the model we found that it is not necessary to restrict the 
parameters b w , 6 2 o an d 6/ (which influence the behaviour of F 2 at x — > 1) to positive values. 



Indeed, there are two domains when x ? 
the x-behaviour of F 2 is controlled by hi 



1. In the first one, Q is large, Q 3> W — m 2 , 
rather than by fojo- in the second domain W 2 is 



small, W 2 



m 



0, it is a region where the Regge approach cannot be applied. 



Thus, we allowed the negative values for b 10 , 6 2 o and bf and found out a new set of 
parameters which, though noticeably different from the old one, leads to a very good x 2 ~ 
1.07. 



3 Extracting the local slopes B x (x, Q 2 ) from the SF 
data. 

As noted in the introduction, the x— slope is a precious tool to settle if, either yes or no, a 
damping of the HERA effect does exist. Actually, we need a set of experimental data on 
B x (x,Q 2 ), at fixed and not too high x versus sufficiently high Q 2 . Unfortunately, because 
the x— slope is not a measurable observable it should be extracted, when possible, from the 
available data on the SF. 

For each given Q 2 , we have an insufficient number of x values for which SF data are 
available, to perform an analysis based on independent x— bins and to extract x-slopes 
with a good accuracy. We adapt the so-called method of "overlapping bins" p2| , previously 



intended for analyzing the local nuclear slope of the first diffraction cone in pp and pp elastic 
scattering. 

Provided that the SF has been measured for a given Q 2 at N x-points lying in some interval 
[x min , x max ], we adopt the following procedure. First, we divide this interval into subintervals 
or elementary "bins" (with rib measurements in each of them, assumed for simplicity to be 
the same for all bins). When the first bin is chosen, the second bin is obtained from the first 
by shifting only one point of measurement (of course one could shift bins by any number of 
points less or equal rib, the shift of one point is the minimal one giving rise to the maximal 
number of overlapping bins). The third bin is obtained from the second one by the shift of 
one point etc... Thus, we define N — rib + 1 overlapping bins for a given Q 2 . For each (fc-th) 
bin, rib must be large enough and its width (in x) small enough to allow fitting the SF with 
the simplest form directly involving the x— slope 

F 2 (x) = A[-\ , (for a given Q 2 ) . (6) 

The parameter B represents the value of the x— slope B (< x >k,Q 2 ), "measured", at Q 2 
and at the "weighted average" < x >k defined in the k-th bin as (see p4j ) 

<x> k =exp{-—^), ke[l,N-n b + l], (7) 

where X{ is the value of x at which the structure function yi is measured with the uncertainty 
Ayf, the summations run over all data points, i = 1,2,..., rib of the bin. This yield the 
"experimental" values of B x (x,Q 2 ) with the corresponding standard errors determined in 
the fit of (6) to the data. Then the procedure is to be repeated for all bins and ultimately 
for the other Q 2 's at which the SF have been measured. 

The next step in extracting and analyzing B x (x, Q 2 ) is the determination of the slopes at 
fixed x as function of Q 2 , making use of results of the first step. As a rule, the sets of < x >&, 
at different Q 2 , do not coincide. So in order to get the x— slope at fixed x and at different Q 2 
we interpolate (or extrapolate but not far from the rc-interval under consideration) already 
extracted B x at the given Q 2 to the chosen x. It can be made assuming a linear x— dependence 
of B x . 



The above mentioned method of overlapping bins is applied to the whole available data 
set. The HERA data (of ZEUS and HI) are analyzed all together [18-25] but independently 
(separately) from older data of the other experimental Collaborations [26-30] because there 
are large gaps between the x-values for these two groups of data (remind that our aim is to 
extract the local B x rather than averaged ones). 

The resulting values of B x (x,Q 2 ) are shown in the Figs .2 (a-e), where they are com- 
pared to the predictions of theoretical models. The results of the interpolation for x = 
0.005, 0.01, 0.05 and 0.08 are presented in the Fig. 3. 

We would like to comment some "technical" points in our analysis and results. 

i) In order to keep its local character to the x-slope and to obtain a maximal possible 
number of "measurements" , we have considered only the cases with four and five points 
in each elementary bin. One can see from the presented figures a quite weak dependence 
of the results on a change in the width of elementary bin. 

ii) In spite of a high accuracy of the recent data from HERA, the dispersion of the SF- 
values strongly influences the resulting values of B x and of its errors. For some bins, 
for example, we were unable to fit (6) to the data with x 2 < 1, so we could not obtain 
in all cases reasonable errors. Nevertheless we show these extracted values (only with 
X 2 < 3)in the figures and use them for interpolation to the x under interest because 
even with these points the extracted set of data for B x (x, Q 2 ) at fixed x is quite poor. 
Of course, this reduces the reliability of our results, but only slightly because the 
number of "bad" bins remains small. 

iii) "Experimental" values of B x at fixed x shown in Fig. 3 are obtained by a linear inter- 
polation within the two subsets of local x-slopes, extracted from the HERA and from 
the fixed target measurements of SF. 

iv) One can see in Fig. 3 that several points deviate strongly from the groups constituted 
by the other points. This is due to the strong influence of the points (of F 2 as well as 
of B x ) which are at the ends of the x-bins and which also "fall out of a common line" 
(see also item ii) above). To solve this problem, a possible way would be to exclude 
some of them from the analysis; another way would be to enlarge the number of points 
in an elementary bin losing, however, the local character of the extracted x— slope. A 
more detailed analysis of the available data related to more numerous measurements 
of the primary observable, F 2 (x, Q 2 ) would be necessary to obtain a better set of data 
ioiB x (x,Q 2 ). 
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Fig. 2 (a-e). The x— slopes B x (x,Q 2 ) at 
fixed Q 2 and versus x, extracted from the 
experimental data on the structure func- 
tion by the method of overlapping bins. 
Triangles (n& = 5) and circles (rib = 4) cor- 
respond to the slopes extracted from the 
HERA data. Stars (n b = 5) and diamonds 
(rib — 4) are the slopes extracted from 
fixed target data. Curves are predictions 
in the DP (solid line), ALLM (dashed line) 
and LKP (dashed-dotted line) models. 
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Fig. 3. The slopes B x interpolated to some x-s as function of Q 2 . The vertical dotted lines 
show the HERA kinematical limits for given x. See Fig. 2. 

4 Results and discussion 

The calculated x— slopes within the three models are in good agreement with the values 
extracted from the experimental SF via the present analysis (in the data bars, see Figs .2), 
except at largest x where the predictions are slightly higher than experimental points. 

Except for the highest Q 2 , in Fig. 2(e), and for the lowest x in all figures, the theoretical 
curves almost coincide. To that respect, the slope does not appear as a good criterion for 
discriminating among the various models. However, we found that the three models predict 



10 



noticeably different x— slopes when x is very small, outside reported experimental range (the 
difference growing when x decreases) Q. For example, we find B x (x = 10~ 4 , Q 2 = 1 OOGeV 2 ) ~ 
0.2 for DP, « 0.4 for LKP, and « 0.3 for ALLM. 

Thus, an extension of the x, Q 2 range for B x towards small x, requiring much effort in 
measuring and analyzing the FS, would be very interesting to seek the experimental evidence 
of the validity (or non validity) of the Froissart-Martin bound in deep inelastic scattering by 
ensuring whether (or not) B x {x —> 0,Q 2 ) — > 0. Furthermore, exploration of this kinematical 
region allows to select a model upon experimental grounds. 

The damping effect is characterized by a slowing down of the x— slope growth, when Q 2 
increases for any x, resulting in a turn-over in the behavior of th x— slope. 

From the theoretical point of view, we recall that we tested only three models || |TB|, |T7| ], 
which assume quite different asymptotic behaviors of the SF and consequently of B x (x <C 



1,Q 2 — > oo) (for more details see [17|). The Fig. 3 shows the x— slope calculated in these 
models and plotted versus Q 2 , for some x— values. Only LKP and DP models do predict [] 
a bump in the curve for Q 2 above 100 GeV 2 , ALLM does not : consequently, we cannot 
rely on theoretical predictions alone due to that model dependence. 

From the experimental point of view, the figure shows data points resulting from our 
analysis, seeming to correspond to a maximum. However, the most interesting part of the 
Q 2 range to explore in the future (because situated in the downhill of the predicted x— slope 
and allowing to conclude unambiguously) lies above the actual measurements and not too 
far from the actual kinematical limit (at least for HERA facilities) . 

We conclude that a deeper theoretical knowledge of the high energy positron-proton in- 
teraction or/and more precise measurements and analysis, opened in higher Q 2 range, are 
needed to confirm or to disprove without ambiguity a bump structure in the behavior of 
B x (x,Q 2 ) at small fixed x. 
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